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Abstract 

Two experiments have been conducted to acquire 
data for the validation of computational fluid dynamics 
(CFD) codes used in the design of supersonic combustors. 
The first experiment is a study of a supersonic coaxial jet 
into stagnant air in which the center jet is of a light gas, 
the coflow jet is of air, and the mixing layer between them 
is compressible. The jet flow field is characterized using 
schlieren imaging, surveys with Pitot, total temperature 
and gas sampling probes, and RELIEF velocimetry. 
VULCAN, a structured grid CFD code, is used to solve 
for the nozzle and jet flow. The second experiment is a 
study of a supersonic combustor consisting of a diverging 
duct with single downstream-angled wall injector. 
Entrance Mach number is 2 and enthalpy is nominally that 
of Mach 7 flight. Coherent anti-Stokes Raman 
spectroscopy (CARS) has been used to obtain nitrogen 
temperature in planes of the flow, and surface pressures 
and temperatures have also been acquired. Modern- 
design-of-experiment techniques have been used to 
maximize the quality of the data set. 

Nomenclature 

k s Thermal conductivity of wall 

n s Number of samples 

n p Number of parameters 

Pamb Ambient pressure 

Pexit Nozzle exit pressure 

P ref cj Center-jet nozzle reference pressure 

P ref cofiow Coflow nozzle reference pressure 

Pr T Turbulent Prandtl number 

q Surface heat flux 

Sc t Turbulent Schmidt number 

t Time 

T Temperature 

T amh Ambient temperature 

T t C j Center-jet nozzle total temperature 

Tf, cofiow Coflow nozzle total temperature 


u Velocity 

x, y, z Position coordinates 
a s Thermal diffusivity of wall 
X Mole fraction center-jet gas 
G Standard deviation 

Introduction 

Computational fluid dynamics (CFD) codes are 
extensively employed in the design of high-speed air. 
breathing engines. CFD analysis based on the Reynolds 
averaged Navier-Stokes equations uses models for the 
turbulent fluxes that employ many ad hoc assumptions and 
empirically determined coefficients. Typically, these 
models cannot be applied with confidence to a class of 
flow for which they have not been validated. Two studies 
have been conducted to provide data suitable for code 
development and testing. 

The first experiment 1 ' 2 ' 3 is a study of a coaxial jet 
discharging into stagnant laboratory air, with center jet of 
a mixture of 5% oxygen and 95% helium by volume and 
coflow jet of air. The exit flow pressure of both center-jet 
and coflow nozzles is 1 atmosphere. The presence of 
oxygen in the center jet is to allow the use of an oxygen 
flow-tagging technique (RELIEF 4 ) to obtain non-intrusive 
velocity measurements. Both jets are nominally Mach 1 .8, 
but, because of its greater speed of sound, the center jet 
flow has more than twice the velocity of the coflow. The 
mixing layer which forms between the center jet and the 
coflow near the nozzle exit is compressible, with a 
calculated convective Mach number 5 of - 0.7. 

This geometry has several advantages: The 
streamwise development of the flow is generally 
dominated by turbulent stresses (rather than pressure 
forces), and thus calculations are sensitive to turbulence 
modeling. It includes features present in supersonic 
combustors, such as a compressible mixing layer near the 
nozzle exit and a light-gas/air plume downstream. Since it 
is a free jet, it provides easy access for both optical 
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instrumentation and probes. Since it is axisymmetric, it 
requires fewer experimental measurements to fully 
characterize, and calculations can be performed with more 
modest computer resources. However, weak shock waves 
formed at the nozzle exit strengthen and turn normal as 
they approach the axis, complicating the flow. Care is thus 
taken in the design of the facility to provide as near as 
possible to 1 -D flow at the exit of both center and coflow 
nozzles, and to minimize the strength of waves generated 
at the nozzle exit. 

Results from this experiment are compared to CFD 
solutions obtained by VULCAN, a previously developed 
code used in engine analysis. 6 

The second experiment 7 is a study of a supersonic 
combustor consisting of a diverging duct with single 
downstream-angled wall injector. Thus, the geometry is 
relatively simple and large regions of subsonic 
recirculating flow are avoided. The nominal entrance 
Mach number is 2 and the enthalpy of the test gas (hot air 
“simulant”) is nominally that of Mach 7 flight. It was 
believed, on the basis of calculations performed 8 that this 
would produce mixing-limited flow, that is to say, one for 
which chemical reaction to equilibrium proceeds at a 
much greater rate than mixing. It 
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Application of a complicated technique like CARS 
in a high-speed engine environment is not routine. Since it 
is a pointwise (rather than planar) technique, building a 
“picture” of the internal temperature field of the 
combustor requires hundreds of facility runs, which is 
expensive. Thus, modern-design-of-experiments (MDOE) 
techniques are used to minimize the quantity of data 
required to meet the goals of this work and to minimize 
systematic errors associated with random errors. Details 
of the MDOE aspects are not discussed in this paper, but 
may be found in Ref. 7. 



The coaxial jet assembly is shown in Fig. 1 . It is 
axisymmetric and consists of outer and center bodies. The 
passage formed by the space between these bodies, and by 
the interior passage of the center body, are nozzles 
designed by the method of characteristics to produce 1 -D 
flow at their exits. 

The nozzle assembly is joined to the Transverse Jet 
Facility, located in the laboratories of the Hypersonic 
Airbreathing Propulsion Branch at NASA Langley 


later proved that this was not the 
case. 

The primary experimental 
technique employed is coherent 
anti-Stokes Raman spectroscopy, 
known by its acronym CARS. An 
introduction to CARS is given by 
Eckbreth 9 , and an application of 
CARS to supersonic combustors is 
given by Smith et al. 10 The species 
probed is molecular nitrogen and 
the quantity measured is 
temperature. Intrusive probes, such 
as Pitot, total temperature, hot- 
wire, etc., are not used due to 
access difficulty and high heat flux 
in the combustor, and because they 
may alter the flow. CARS has 
several advantages over other 
optical methods. It is a relatively 
mature and well-understood 
technique. Signal levels are 
relatively high and the signal is in 
the form of a coherent (laser) beam 
that can be collected through small 
windows. Incoherent (non-CARS) 



interferences are rejected by spatial 


All dimensions in mm 


filtering. 


Figure 1. Coaxial jet flow facility. 
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Research Center. The plenum of this facility contains 
porous plates for acoustic dampening and screens for flow 
conditioning. Air is provided to the facility from a central 
air station, and the helium-oxygen mixture is provided to 
the center body from a bottle trailer containing premixed 
gas. 

The assembly is instrumented with pressure taps: one 
in the center body just downstream of the screens, one in 
the facility plenum, and one in the outer body near the exit 
of the coflow nozzle (in a region where the flow has 
reached its exit condition). Thermocouples are located in 
the gas supply lines. Ambient (barometric) pressure and 
ambient temperature are read. The values of these various 
quantities during the probe surveys, and their respective 
uncertainties (95% probability band) are given in Tab. 1. 
Uncertainties are due to facility unsteadiness and 
variations in set point, and do not include ±0.5% in 
pressures and ±2 K in temperatures due to transducer 
error. 


P ref coflow (kPa) 

580 ±2 

T t ,coflow (K) 

300 ±6 

Pref Cj/p ref coflow 

1.060 ±0.008 

T coflow 

1.02 ±0.05 

P amt/p ref coflow 

0.1758 ±0.0012 

Tamt/T ^ coflow 

0.982 ±0.017 

Pexi/P ref coflow 

0.1748 ±0.0005 


Table 1. Experimental flow parameters. 

Flow Field Measurements 

Various types of flow field measurement have been 
performed. The flow has been visualized with 
conventional schlieren and shadowgraph. Pitot, gas 
sampling, and total temperature probes have been 
employed to survey the flow. Probe survey locations are 
listed in Tab. 2, and are also shown in Fig. 3(b). 
References 1 and 2 give details of these measurements. 

Survey probe tips are cylindrical and cut square, with 
outside/inside diameters respectively of the Pitot probe 
0.64 mm/0.36 mm, and of both gas sampling probe and 
total temperature probe 1.27 mm/0.76 mm. The gas 
sampling probe and tubing internal diameters are sized to 
avoid choking the sample gas flow, ensuring shock 
attachment at the probe tip. The total temperature probe is 
a miniature shrouded, vented thermocouple. The probe 
incorporates a commercial microminiature thermocouple 
junction at the tip of a 0.20 mm diameter “needle”. Errors 
in Pitot pressure due to pressure transducer error are 
±0.5%. Error in total temperature due to thermocouple 
error is ±2 K. In addition, the total temperature probe was 
found to read about 1% low, due to incomplete stagnation 
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of the flow at the sensor and/or radiation losses. 


Number 

x, mm 

x, mm 
RELIEF 

1 

0.13 

2 

2 

3.1 

5 

3 

10.0 

12 

4 

17.9 

17 

5 

27.8 

27 

6 

42.9 

42 

7 

61.9 

62 

8 

81.1 

82 

9 

100.6 

102 

10 

121.4 

123 

11 

150.8 

153 

12 

181.0 

190 

13 

220.4 

220 

14 

261.0 

258 


Table 2. Experimental survey locations. 

The mole fraction of the center-jet gas (i.e., the 
helium-oxygen mixture) in the gas withdrawn from the 
flow, /, is found in real time by a hot-film probe based 
system 11 . The largest contribution to the uncertainty of the 
system is the manufacturer-quoted ±1 % of full scale in the 
mass flow controller used to provide a known helium- 
oxygen/air mixture to calibrate the system. Maximum 
uncertainty in mole fraction of helium-oxygen is in the 
range ±1-1.5%, but uncertainty is less than this for mole 
fractions around 0.0 or 1.0, where there is no uncertainty 
in the composition of the calibration mixture. 

The probes were mounted in a diamond-airfoil strut, 
and translated in the flow by a two-component stepping- 
motor driven translation stage. Probe “zero” location was 
determined using machined fixtures mounted to the nozzle 
exit (conical extension cap removed). Surveys were 
conducted across a diameter of the flow. Analysis of the 
data to find the best-fit center showed it to be within 0.4 
mm (95% of the time) of the measured center. Thus, 
probe surveys are taken to pass through the axis of the jet 
±0.4 mm. Survey data presented have been shifted (by 
less than ±0.4 mm) so that the best fit center lies at y = 0. 
Resulting data are found to be almost perfectly 
symmetrical. 

In addition to these “conventional” techniques, the 
RELIEF 3 4 (Raman Excitation plus Laser-Induced 
Electronic Fluorescence) oxygen flow tagging technique, 
illustrated in Fig. 2, has been used to provide 
measurements of (instantaneous) axial component 
velocity. RELIEF is a time-of-flight technique that 
involves two steps. In the first (or “tagging”) step, oxygen 
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in a line segment of the flow is 
excited to a non-equilibrium 
vibrational state by stimulated 
Raman scattering. This is 
achieved by focusing (50 cm 
focal length) collinear laser 
beams at 532nm and 580 nm. 
These beams are generated by 
passing a 200 mJ doubled 
Nd:YAG laser beam (532 nm) 
through a 6.9 MPa Raman cell 
containing a 50:50 mixture of 
helium and oxygen. The Raman 
cell is seeded with light from a 
broadband dye laser pumped by 
doubled residual infrared light 
from the Nd:YAG laser. The 
non-equilibrium oxygen returns 
to equilibrium only slowly as it 
convects with the flow. In the 



second "probe” step of the 
technique, the non-equilibrium region is found by laser- 
induced fluorescence imaging. This is achieved with a 20 
mJ narrow band (approximately 0.5 cm 1 ) ArF excimer 
(193 nm) laser beam cylindrically focused to a 10 mm 
high x 0.5 mm thick sheet in the region where the tagged 
flow is expected to be. The resulting fluorescence is 
imaged using a double intensified video-rate CCD camera, 
with f/4.5 UV lens and extension rings for closeup 
operation. Data were acquired at 5 Hz. 

The resulting data consist of images of displaced 
line pairs, acquired either at different delay times after the 
tag, or with one of the lines acquired prior to operation of 
the jet (i.e., with zero flow velocity). The instantaneous 
velocity is determined by finding, in subsequent data 
reduction, the line displacement at various points along it, 
and dividing by the probe delay time. A calibration is 
required to establish the relationship between position in 
the image and position in space. Mean w-component 
velocity and root mean square fluctuation have been 
obtained by this technique. Uncertainties in this data are 
approximately ±3% due to uncertainty in the 
magnification factor between flowfield and image, and 
uncertainty in the zero point. 


distributed among five blocks. Grid points are clustered 
near the walls of the nozzles to resolve the boundary 
layers, at the exit of the center-jet nozzle to resolve the 
recirculation zone and shocks in the vicinity of the nozzle 
lip, and to a lesser degree near the axis to resolve shock 
reflections. The distance from the wall of the centers of 
the closest cells is less than y + =1.5 for all surfaces. 

The walls are specified to be adiabatic and no slip. 
Total pressure and temperature conditions are specified at 
subsonic inflow/outflow planes, while the code switches 
to extrapolation where the code detects that outflow is 
supersonic. The flow is assumed to be axisymmetric. At 
the exterior boundary the composition is air with density 
of 1.177 kg/m 3 and pressure (p amb ) 101.3 kPa. At the 
coflow nozzle inflow boundary the composition is air with 
total density 6.735 kg/m 3 and total pressure ( p re f C ofioJ 
580.0 kPa. At the center-jet nozzle inflow boundary the 
composition is 0.7039 by mass He and 0.2961 by mass 0 2 
with total density 1.3343 kg/m 3 and total pressure 628.3 
kPa (computed from p refC j and the area ratio between the 
reference plane and sonic throat, assuming quasi- 1-D 
flow). 

The flow is assumed to be turbulent, and Wilcox’s 12 


Calculation 

The Favre-averaged Navier-Stokes equations are 
solved using VULCAN, a structured, finite-volume CFD 
code. The calculation assumes an axisymmetric flow of a 
mixture of thermally perfect gases. The calculation was 
performed on a structured grid generated by a separate, 
commercial code. There are a total of 188,080 cells, 


k-cb turbulence model is used with the high Reynolds 
number model. The compressibility correction proposed 
by Wilcox was not used, but Wilcox’s generalization of 
Pope’s modification to the k -e model (which attempts 
to resolve the "round jet/ plane jet anomaly”) was. 
Turbulent Prandtl number and Schmidt number were set 
equal to 0.75. More details of the calculation may be 
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Figure 3. (a) (left) Schlieren image with vertical knife edge (conical 
extension cap removed) and (b) (right) computed Mach number 
distribution with data survey lines. 


found in Ref. 3. 

Results 

Figure 3(a) is a typical schlieren 
image (with knife edge vertical) 
showing the jet with coflow nozzle 
conical extension ring removed. 

Vertical dark and bright bands are 
due to transverse gradients of 
refractive index. Notice the shock- 
expansion wave structure emanating 
outward from the (0.25 mm thick) 
center-body lip. Similar waves 
propagate in the center jet, but are 
not visible in the schlieren due to the 
low refractive index there. The 
continuation of these initially inward 
propagating waves, after they have 
crossed at the axis and passed out of 
the center jet into the coflow air, is 
visible. 

Figure 3(b) is a flooded contour 
plot of Mach number from the CFD 
calculation. The results may be 
qualitatively compared to the 
schlieren. The waves seen radiating 
from the center-jet nozzle lip in the 
schlieren are found in the 
calculation, though are not fully 
resolved. A more detailed inspection 
shows that the wave from the center- 
jet nozzle forms a normal shock 
where it intersects the axis. This 
results in a small deficit in total 
pressure at the axis that is visible 
downstream of the shock in both 
CFD and experimental Pitot pressure (see Fig. 6). This 
deficit persists as far downstream as Plane 9 before being 
obscured by the mixing of the coflow into the center jet. 

Figures 4 - 8 show comparisons between the results 
of the experiment and the results of the CFD calculations. 
The range of y in the plots does not correspond to the full 
range of the data or of the calculation, but is truncated to 
show more clearly the regions of interest. In these figures, 
y is given in m and u in m/s. 

The mole fraction centerjet gas data is shown in Fig. 
4. The centerjet spreads smoothly, with the peak^falling 
below 1.0 downstream of Plane 11. The experimental 
values are well reproduced by the calculation near the 
axis, but, moving away from the axis, the calculation is 
first high and then, near^ = 0, too low. The calculation is 


discontinuous in slope at x = 0 (a most unphysical 
behavior). This discontinuity cannot be attributed to 
inadequate grid resolution or extent of domain, and is 
believed to be a problem with the turbulence model. 

The mean velocity data is shown in Fig. 5. At Plane 1 
there is a layer with velocity deficit at the boundary 
between the centerjet and the coflow that is several times 
the thickness of the nozzle lip. This layer results from the 
merging of the coflow nozzle inner surface boundary layer 
and the region of separation at the lip. The spikes in 
velocity near the edge of the centerjet are due to the shock 
waves emanating from the nozzle lip. Downstream of the 
nozzle exit the velocity data is consistent with the^data, 
showing a similar spread of the centerjet, but the drop in 
peak velocity below the nozzle exit value begins 
significantly further upstream then the drop in ^ nearer 
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Plane 7. Calculated velocity is high compared to the data 
near the axis, but near the edge of the centerjet it is low 
or, further downstream, close to the data. 



y 


Figure 4. Mole fraction center-jet gas at several data 
planes. 



Figure 5. Mean velocity. 

The Pitot pressure is shown in Fig. 6. At Plane 1, as 
with velocity, there is a layer with reduced Pitot pressure 
at the boundary between the centerjet and the coflow. 
Small axisymmetric irregularities in Pitot pressure in the 
centerjet (-0.005 m < y < 0.005 m) may be attributed to 
machining flaws in the center-jet nozzle. In general, 
however, experiment and calculation agree very well, 
indicating that the calculation of the flow in the nozzles 
was good. Downstream of Plane 1 the centerjet spreads, 
with Pitot pressure near the axis falling downstream of 
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Plane 10 and rising in the wake of the nozzle lip. As with 
X, erroneous discontinuities in slope may be observed in 
the calculation at the outer boundary of the centerjet and 
the inner boundary of the outer (coflow-ambient) shear 
layer. 



y 


Figure 6. Pitot pressure. 

Total temperature was acquired at Plane 9 only 
(Fig. 7). On the axis and in the coflow experimental 
data are about 1 % below the known supply gas 
temperatures, due to probe error. Otherwise, the 
calculation agrees well with the experiment, 
reproducing both overshoot and undershoot in the data. 



y 


Figure 7. Total temperature at Plane 9. 
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Supersonic Combustor Experiment 

Flow Facility 

This experiment was conducted in NASA Langley’s 
Direct-Connect Supersonic Combustion Test Facility 
(DCSCTF) 13 . “Vitiated air” is produced at high pressure 
in the “heater”, shown in Fig. 8. Oxygen and air are 
premixed and hydrogen is burned in the mixture. Flow 
rates are selected so that the mass fraction of oxygen in 
the resulting products is the same as that of standard air. 
The test enthalpy is nominally that of Mach 7 flight. The 
vitiated air is accelerated through a water-cooled 
convergent-divergent nozzle and enters the test model. 

Gas flow rates to the heater are: 0.915±0.008 kg/s 
air, 0.0284±0.0006 kg/s hydrogen, and 0.300±0.005 kg/s 
oxygen. The heater stagnation pressure is 0.765±0.008 
MPa. These uncertainties are due to the random run-to- 
run variations and do not include a ±3% uncertainty in the 
mass flow rate measurements. 

Heater and nozzle exit conditions are estimated from 
the flow rates, heater pressure, nozzle minimum area and 
exit area using one-dimensional (ID) analysis 14 . The flow 
exiting the heater into the nozzle is assumed to be in 
thermodynamic equilibrium, but has unknown enthalpy 
due to heat lost to the structure and cooling water. 
Enthalpy is found from the known mass flow rate, 
geometrical area of the nozzle (sonic) throat, assuming 
isentropic flow in the nozzle and 1-D flow at the throat. 
Nozzle exit conditions are computed similarly from the 
geometrical exit area. Calculations assuming equilibrium 
and frozen composition in the nozzle differ in minor 
species concentration, but not significantly in major 
species, temperature or pressure. The nominal calculated 
conditions, and uncertainties due to mass flow rate 
measurement error and run-to-run variations in heater 
conditions are: heater stagnation temperature 1827+75 K, 
exit temperature 1187+60 K, exit pressure 100+1.5 kPa, 
exit Mach number 1.989+0.005. Errors arising in the 
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calculation due to the assumption of ID flow (the 
effects of non-uniform composition, boundary layers, 
etc.) are not considered. 

A study of the flow at the exit of the facility 
nozzle was conducted previously 15 . A Pitot probe 
rake was employed to map the exit Pitot pressure. 
Additionally the flowfield at the exit of the nozzle 
was visualized. Silane (SiH 4 ) was added to the heater 
hydrogen and burned to form silica particles in the 
heater. The particles were illuminated by a pulsed 
laser-sheet and imaged with a CCD camera. Results 
were compared to CFD calculations of the nozzle 
flow. The flow at the nozzle exit was not completely 
ID, but the computed Pitot pressure distribution 
agreed well with measurement. The flow appeared 
well mixed. 

The test model is shown in Fig. 9. There are two 
main sections: the copper section upstream and the carbon 
steel section downstream. Stainless steel flanges and 
carbon gaskets separate these sections from each other 
and the nozzle. The internal passage, from left to right, 
has a constant area segment, a small outward step at the 
top wall, a second short constant area segment followed 
by a constant 3° divergence of the top wall. The span is 
constant at 87.88 mm. Five small pilot fuel injector holes 
are located ahead of the step, and the main fuel injector is 
located just downstream of the start of the 3° divergence. 
The injection angle is 30° to the opposite wall. The 
injector nozzle is designed by the method of 
characteristics to produce Mach 2.5, ID flow at the 
injector exit. Hydrogen injection is provided at a pressure 
of 2.12+0.07 MPa, temperature of 302+4 K, and 
equivalence ratio of 0.99+0.04. On some runs, additional 
hydrogen injection is provided by the 5 pilot injectors at 
the same nominal temperature and a total equivalence 
ratio of 0. 148 ±0.008. The pilots are turned on and off at 
the same time as the main fuel injector. 

The duct is uncooled; however, the wall thickness of 
the copper duct is greater than 32 mm and the carbon steel 
duct wall thickness is 19 mm. Thus, run times fueled in 
excess of 20 s are possible (and much greater if unfueled). 
With atmospheric temperature air flowing in the model 
between runs, runs could be repeated every 10 - 15 
minutes. 

The model is equipped with 7 slots to allow the 
CARS beams to penetrate the duct, of which slots 1, 3, 5, 
6, 7, depicted in Fig. 9(a), are used in this study. The slots 
are in pairs, one on each side of the duct, 4.8 mm wide, 
extending the full height of the duct. When not in use the 
slots are plugged flush to the wall. Windows covering the 
slots are mounted at the end of short rectangular tubes. 
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(a) 




Figure 9. Test model: (a) nozzle, copper and carbon steel (C.S.) duct 
sections, (b) detail in vicinity of fuel injector and pilots. 


The model is additionally instrumented with both 
pressure taps and wall temperature probes. More details of 
this instrumentation may be found in Ref. 7. 

CARS Technique 

The CARS system uses an unseeded Spectra-Physics 
DCR-4 pulsed Nd: YAG laser, frequency doubled to 532 
nm. The nominal power at 532 nm is 550 mJ per pulse 
and repetition rate is 10 Hz. A broadband dye laser 
utilizing two longitudinally pumped Brewster’s angle 
pumped dye cells is employed in the system. Dye laser 
wavelength is centered between 605 nm and 606 nm to 
match the Raman shift of nitrogen by adjusting the dye 
concentration. The dye laser and two 532 nm beams are 
combined at a dichroic mirror and relayed via a periscope 
to a spherical lens. The three beams are crossed at their 
focal points in a vertical planar BOXCARS 
configuration 3 . 

At the focus, the diameters of the 532 nm and dye 
beams are respectively - 0.12 mm and - 0.15 mm. The 
length of the measurement volume is found by translating 
the CARS measurement volume through a thin planar jet 
of nitrogen surrounded by a coflowing jet of helium. The 
length over which CARS signal is recorded is ~ 4.5 mm 
and the full width half maximum (FWHM) of the signal 
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distribution is ~ 2.25 mm. The beam 
energy levels per pulse obtained at the 
focusing lens are -85 mJ for each green, 
and from 12 mJ to 24 mJ for the dye. 

The beams (including the CARS 
signal beam) are relayed via a second 
spherical (collimating) lens and a second 
periscope back to the optical bench. The 
CARS beam is separated, directed 
through additional filters as needed and 
a polarizer that allows only horizontally 
polarized light to pass, then focused to 
the entrance of a 1 m monochrometer 
with 1200 groove/mm grating. An 
EG&G PAR model 1420 intensified, 
linear, self-scanned silicon photodiode 
array detector (IPDA) is mounted at the 
exit plane of the detector. The detector 
consists of 1024 elements of which the 
central 598 elements are used. An 
optical splitter 16 creates a secondary 
signal on the detector, identical to the 
primary but offset by 290 pixels and 
6. 1 % the intensity. When the intensity of 
the primary signal exceeds the dynamic 
range of the detector, the secondary 
signal is used for analysis. 

The two top prisms of the periscope are mounted on 
stepping motor driven vertical translation stages. The two 
bottom prisms and the vertical translation stages are 
mounted on similar horizontal stages. By translating the 
vertical and/or horizontal stages in tandem the 
measurement volume could be moved in the y and/or z 
direction. CARS data acquisition is under the control of a 
personal computer (PC). Two types of acquisition are 
employed. In the first, data is acquired at a single point in 
space. In the second, data is acquired while either the 
vertical or the horizontal stages are in constant velocity 
motion. 

CARS data are acquired in the supersonic combustor 
during multiple sets of test runs. During a set of runs 
(which might last as long as 5 hours), access to the model 
and optical system is prohibited for safety reasons. Test 
runs consisted of approximately 5 s during which the 
heater is operating but no fuel is injected in the model, 
followed by from 11 s to 20 s during which fuel is 
injected. CARS data is acquired over a period 2 s shorter 
than the period of fuel injection. Immediately after a run, 
10 s of data is acquired with the system operating as 
before, the dye laser beam blocked by a remotely operated 
beam block. These “background” scans measure non- 
CARS interferences such as scattered laser light. 
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CARS data are analyzed on a separate workstation. 
Prescans are subtracted from data scans. Background 
scans (after subtraction of prescans) are averaged and 
subtracted from data scans. Both primary and secondary 
(produced by the splitter) CARS signals are contained 
within the data scan. If the primary is saturated, the 
secondary is selected for analysis. Data scans are divided, 
pixel by pixel, by the reference spectrum to remove the 
effect of the dye laser spectral power distribution, and 
normalized to unit area (primary or secondary). Data are 
compared to a library of similarly normalized theoretical 
spectra to determine the temperature and nitrogen 
concentration. The pixel location of the start of the 
theoretical spectra is allowed to vary for best fit. The 
combination of temperature, concentration, and pixel 
location that produces the least mean square deviation 
between theory and data is selected. 

Theoretical CARS spectra are generated using the 
program CARSFT 17 . The combustion gases are assumed 
to be a mixture of nitrogen and non-resonant buffer gas, 
both having non-resonant susceptibility of 8.5xl0~ 18 
cm 3 /erg. The static pressure is assumed to be 1 
atmosphere, although, in reality, the pressure varied (see 
Fig. 10). The Exponential Gap Model for collisional 
narrowing of the Raman line shape is used. A 532 nm 
laser line width of 1 cm 1 is assumed. An experimentally 
determined instrument probe function is used. The 
reference CARS spectrum can shift significantly in center 
wavelength during a set of runs. Techniques developed to 
derive suitable reference spectra from the CARS data are 
described in Ref. 7. 

The techniques used for acquisition and analysis of 
CARS data in the supersonic combustor were tested in a 
“Hencken” adiabatic, flat-flame burner burning hydrogen 
in air. Equivalence ratio (ratio of hydrogen rate to 
stoichiometric hydrogen rate for given air flow) was 
varied and the measured temperature compared to 
calculations based on measured flow rates and equilibrium 
chemistry (including minor species). Tests were 
conducted in which the total laser power was varied from 
200 mJ to 550 mJ. Also data were obtained in which, 
through the use of different neutral density filters, the 
signal in the primary is saturated, forcing use of the 
secondary in analysis. No trends are found with either of 
these variables, indicating that the nitrogen spectrum is 
not saturated by high laser powers and that the splitter 
device works well. The average of all the data at an 
equivalence ratio 1.0 is 2360 K, compared to the 
theoretical value of 2380 K. The measurements agreed 
with calculation within ±100 K at an equivalence ratio 
less than or equal to one. 


Surface Pressure and Temperature 

Surface pressure and temperature data are presented 
for two typical runs, one in which the pilot injectors are 
operating and one in which they are not. These runs are 
ones for which the gas flow rates to the heater, injector 
and pilot and the heater pressure are all very close to their 
respective averages over the total set of runs. 

Surface pressure distributions for the pressure taps at 
the bottom wall centerline, averaged over 1 s intervals, are 
shown in Fig. 10. The heater is initiated at time t 0 = Is 
and fuel injection commences at time t x = 6.0 s (piloted) 
or 6.4 s (unpiloted). Data are shown 10 s and 22 s into the 
run. Pressures vary widely in the upstream region due to 
the complex nature of the shock wave system created by 
the injectors and step. The pressure for the piloted case is 
higher than the unpiloted case between the pilot (x = 
0.074 m) and about x = 0.7 m, due to combustion of gas 
from the pilot and main fuel injector. In the unpiloted 
case, pressure tends to fall moving downstream due to 
divergence of the duct, until 0.5 m where it rises rapidly, 
peaking at about 0.75 m. Presumably, there is minimal 
combustion upstream of 0.5 m. Downstream of 0.75 mthe 
pressure drops smoothly in both cases but is higher in the 
unpiloted case, despite the greater total injected fuel rate 
in the piloted case. Differences between the two cases 
suggest that significant combustion of the fuel does not 
take place upstream of 0.5 m in the unpiloted case, but 
that the combustion then proceeds to completion by 0.75 
m. It is not believed that fuel and air are fully mixed at 
0.75 m, so that further mixing and combustion occurs 
downstream. 



centerline of bottom wall. 

Comparison between measurements at 10 s and at 22 
s reveals only small differences. There is no suggestion 
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that the combustion delay experienced by the unpiloted 
case is affected by the increase in surface temperature that 
occurs during a run. Surface pressure varies as much as 20 
kPa or more between top, bottom, and sidewall taps (not 
shown) upstream of about 0.65 m, consistent with the 
effects of shock waves. Downstream of this point, there 
are no differences between the walls. 

In the copper section, temperature is typically about 
- 360 K at the start of fuel injection but rises to as high as 
~ 610 K. In the carbon steel section, it typically is ~ 440 
K at the start and as high as - 950 K at the end. These 
variations in surface temperature are least-squares fit to 
the solution for wall temperature of a semi-infinite body at 
initially uniform temperature, subject to steps in surface 
heat flux at heater start and fuel injection start. Fit 
parameters include heat flux with heater only, q (h and heat 
flux rise due to fuel injection, q x . The fit is conducted out 
to t = 11 s. Representative temperature histories and fits 
are shown in Fig. 11. Fits often diverge from the data 
beyond 11s, indicating that heat flux continues to change 

slowly during the run. The material property yfa^ / k s is 

taken to be 36.7 kW K/m 2 s 1/2 for the copper duct and 12 
kW K/rrrs 1/2 for the steel duct. These values have not been 
verified, so this analysis should not be relied upon except 
in a relative sense (i.e., case to case, location to location 
within the copper duct, injection to no injection). 



Figure 11. Wall temperature history, unpiloted, top 
wall: data and fits. 


Heat fluxes are presented in Fig. 12, on both top and 
bottom walls, at the centerline. Heat flux varies from 1.0 
MW/m 2 to 0.3 MW/m 2 without fuel injection, and from 
0.7 MW/m 2 to 1.8 MW/m 2 with injection. Without 
injection, heat flux varies at a given location 10% to 30% 
from case to case, reflecting variation in the initial 
temperature of the wall between runs. With injection, heat 
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flux shows large increases relative to before injection. 
Large increases occur in the piloted case on the top wall, 
downstream of the main injector (x = 0.166 m). In the 
unpiloted case, a smaller heat flux rise occurs upstream of 
0.5 m (where the pressure rise starts), indicating either the 
start of combustion near the wall, or an increase in heat 
transfer coefficient. (An increase in heat transfer 
coefficient could occur under the fuel plume due to a 
streamwise vortex pair.) 

2 r 

Unpiloted, q 0 
B - - Bottom 

- -0 - - Top 

Unpiloted, qo+qi 
■ Bottom 

- -• - -Top 
Piloted, q 0 

B Bottom 

0 Top 

Piloted, qo + qi 
■ Bottom 


X, m 

Figure 12. Surface heat flux distributions: top and 
bottom walls, centerline. 

CARS Temperatures 

Data was acquired over 201 facility runs over 10 test 
days. Except for one day, when laser beams clipped the 
edge of the duct window slots due to thermal expansion 
and movements of the duct, the vast majority of the data 
were found acceptable and analyzed. At each plane, data 
were acquired at 6 or 7 fixed points near the horizontal 
centerline. Data were also acquired during 16 s of 
horizontal motion of the translation stages at 5 mm/s, or 
from 9 s to 18 s of 5 or 6 mm/s vertical motion, during 
which time fuel is continuously injected. 

All CARS temperature data are fit to a cosine series 
bivariate function of order 5 at plane 1 or order 6 at the 
other planes, with the number of fit parameters 
respectively 21 and 28. Commercial software was used 18 . 
2000 to 4000 data points were acquired per plane and the 
standard deviation of the data from the fit at the various 
planes ranges from 196 K to 304 K. Thus, the fitted 
functions represent an estimate of the mean temperature 
distribution with mean uncertainty, given by 

1.98 <Jyjn p / n s 19 , from 36 K to 59 K depending on the 

plane. It is important to point out that the uncertainty in 
the surface fits is lower near the center of the 
measurement plane and higher near the edge. This 
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uncertainty does not include any non-random error 
components. 

Figure 13 contains 3-dimensional cutaway views of 
the duct showing contour plots of the fitted temperature 
functions. Flow is from top left to bottom right. Recall 
that the main fuel injector is on the top wall between 
Planes 1 and 3, and the pilot injectors are in the top wall 
upstream of Plane 1. 

At Plane 1 in the unpiloted case the temperature is 
fairly uniform, between 1030 K and 1250 K. The mean 
temperature for all the data points of this plane and case 
is 1 162 K. This mean compares favorably with the value 
computed assuming 1-D flow from the heater, which is 
1187±60 K. In the piloted case, the temperature drops 
slightly close to the top wall where the (cold) pilot fuel is 
injected. There is no indication of combustion in this 
plane. 

At Plane 3, in the unpiloted case, the injected fuel 
plume is a region of low temperature with temperatures 
as low as - 250 K at the center. There is no evidence of 
combustion of the injected fuel. In the piloted case, there 
is a band of hot combustion products close to the top 
wall. The center of main fuel plume may be seen as a 
cool region with temperature as low as ~ 650 K, which is 
greater than in the unpiloted case, suggesting some 
combustion. 

For both cases, Plane 5 is similar to Plane 3. In the 
unpiloted case temperatures have risen near the center of 
the fuel plume to a minimum of - 550 K. In the piloted 
case the minimum is ~ 1250 K and the height of the 
region of hot combustion products near the top wall is 
greater. 

At Planes 6 and 7 (unpiloted case only), 
temperatures have risen abruptly as compared to Plane 5, 
suggesting nearly complete combustion, i.e., combustion 
that has consumed all available oxygen or fuel at a given 
point. The hot region close to the top wall (temperatures 
as high as ~ 2300 K), and to a lesser extent that near the 
bottom, are probably ones in which the fuel-air ratio is 
nearly stoichiometric. The cooler region near the center 
(as low as - 1500 K) is probably fuel rich. Injected fuel 
may not have penetrated to the sidewalls (~ 1200 K to 
1 300 K). Data were not acquired in the piloted case due to 
lack of facility time. However, similarity in the pressure 
distributions downstream of about x = 0.75 m suggest that 
the temperature maps would be at least qualitatively 
similar, except perhaps close to the top wall. 

The various data obtained provide a consistent 
description of the flow field. In the unpiloted case, no 
significant or only small and intermittent combustion of 
the injected fuel is observed ahead of x = 0.5 m. 
Downstream of = 0.75 m combustion appears nearly 


Piloted 

Unpiloted 



Figure 13. Cutaway views of duct showing contours of mean 
temperature for unpiloted and piloted cases. 

complete. In the piloted case, combustion of the pilot fuel 
appears to take place between jc = 0.122 m (at the step) 
and 0.274 m. There also appears to be significant 
combustion of the main injected fuel by this location. 
These results are not consistent with the CFD calculations 
performed prior to the experimental work 8 , which 
predicted combustion in the vicinity of injection in both 
unpiloted and piloted cases. Consequently, this 
experiment provides a test case for CFD that is more 
challenging than anticipated. Accurate calculation will 
require accurate modeling of the chemical kinetics and 
turbulence-chemistry interactions as well as accurate 
modeling of the turbulent mixing. 

Summary 

Two experiments to acquire data for the validation of 
computational fluid dynamics (CFD) codes used in the 
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design of supersonic combustors have been described. 

The first was a study of a supersonic coaxial jet with 
center jet helium and coflowing jet air. Data acquired 
include schlieren visualization, gas sampling, Pitot probe 
surveys, and RELIEF flow tagging. Calculations utilizing 
a structured finite difference code and Wilcox’s k - cb 
model have been presented. The results of the calculations 
showed non-physical discontinuities in slope of mole 
fraction and Pitot pressure that are attributed to 
inadequacies in the turbulence model. 

The second was a study of a supersonic combustor 
with single downstream-angled hydrogen fuel injector. 
Data include CARS temperature maps, and wall pressures 
and temperatures. Modern design of experiments 
techniques were used to maximize data value. It was 
found that (without pilot injectors) ignition did not occur 
until significantly downstream of injection. Previously 
performed CFD calculations suggested that ignition 
occurred close to injection. This discrepancy is attributed 
to inadequacies in the kinetics model and/or the model for 
turbulence chemistry interaction. 

The 1 st author would like to acknowledge the 
support of the NASA Langley Research Center through 
grant NCC 1-370. 
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